Graphite from the viewpoint of Landau level spectroscopy: 
An effective graphene bilayer and monolayer 
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We describe an infrared transmission study of a thin layer of bulk graphite in magnetic fields up to 
B — 34 T. Two series of absorption lines whose energy scales as \fB and B are present in the spectra 
and identified as contributions of massless holes at the H point and massive electrons in the vicinity 
of the K point, respectively. We find that the optical response of the K point electrons corresponds, 
over a wide range of energy and magnetic field, to a graphene bilayer with an effective inter-layer 
coupling 271 , twice the value for a real graphene bilayer, which reflects the crystal ordering of bulk 
graphite along the c axis. The K point electrons thus behave as massive Dirac fermions with a mass 
enhanced twice in comparison to a true graphene bilayer. 

PACS numbers: 71.70.Di, 76.40.+b, 78.30.-j, 81.05.Uw 
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Recent interest in graphene (HQ, a truly two- 
dimensional system with its simple, but nevertheless, 
for solids, unconventional electronic states, has focused 
attention on the properties of Dirac-like fermions in 
condensed matter physics in general. Whereas, two- 
dimensional massless Dirac fermions [H, 13, [Bj] , character- 
istic of graphene have been widely investigated, far fewer 
experiments have been devoted to massive Dirac fermions 
which are specific to a graphene bilayer , which rep- 
resents a further example of a two-dimensional system 
with a highly unusual band structure Perhaps sur- 
prisingly, Dirac dispersion relations can also be found in 
graphite, a three dimensional, bulk material which con- 
sists of Bernal-stacked weakly coupled graphene layers. 

The standard Slonczewski-Weiss-McClure (SWM) 
model of electronic states in graphite 0, [l(| predicts a 
complex form for the in-plane dispersion relation which 
changes considerably depending upon the value of the 
momentum k z in the direction perpendicular to the lay- 
ers. Intriguingly, the SWM model predicts that in 
the vicinity of the H point (k z = 0.5) the in-plane 
dispersion is linear and thus resembles a Dirac cone. 
Such a dispersion has indeed been found in angle re- 
solved photoemission spectroscopy 
spectroscopy 
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14j, as well as in Landau level (LL)- 
spectroscopy [15| . The latter experiments, mainly fo- 
cused on transitions between LLs whose energy scales as 
\f~B, are generally believed to exhibit far richer spectra 
in comparison to true graphene 3, 17 , 3] , reflecting the 
inherent complexity of the SWM model which includes 
no fewer than seven parameters 1||, 2fj| . 

We show in this Letter that infrared magneto- 
absorption spectra of graphite, measured over a wide 
range of the energy and magnetic field, can be interpreted 
in a very simple, transparent and elegant manner. Our 
results confirm, in agreement with theoretical consider- 
ations 2l|, that graphite can be viewed as an effective 



graphene monolayer and bilayer. This theoretical pic- 
ture is derived using a drastically simplified SWM model, 
which includes only two parameters 70 and 71 , describing 
the intra- and inter-layer tunneling respectively. In this 
simplified picture, the dominant contribution to the op- 
tical response is provided by the H point, where electron 
states closely resemble graphene but with an additional 
double degeneracy, and by the K point, where the en- 
ergy spectrum resembles a graphene bilayer, but with an 
effective coupling of 271, twice enhanced compared to a 
real bilayer system. 

Remarkably, using this simple graphene monolayer 
plus bilayer view of graphite, we are able to correctly 
reproduce the magnetic field evolution of all observed 
inter-LL transitions using only the SWM parameters 70 
and 71, with values which perfectly match those derived 
from studies of real graphene monolayer and bilayer sys- 
tems. Interestingly, the electronic states at K point of 
graphite are found to mimic those of the graphene bi- 
layer, but with a doubled value of the effective mass, so 
that they might be useful to further explore the interest- 
ing physics of massive Dirac fermions. 

Thin samples for the transmission measurements were 
prepared by exfoliation of a natural graphite crystal as 
described in Ref. [HI]. Some data is also presented for 
highly oriented pyrolytic graphite, which shows prac- 
tically identical, although slightly less pronounced fea- 
tures, in the magneto-transmission spectra (22|. All 
experiments were carried out on macroscopic, roughly 
circular-shaped samples, of several millimeters in diame- 
ter. Measurements were performed in the Faraday config- 
uration with the magnetic field applied along the c-axis of 
the graphite. All spectra were taken with non-polarized 
light. To measure the magneto-transmittance of the sam- 
ple in the spectral range of 10-700 meV, the radiation of 
globar was delivered via light-pipe optics to the sample. 
The radiation, detected by a Si bolometer, placed directly 
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FIG. 1: (color online) Transmission spectra of a thin graphite 
layer as a function of the magnetic field in the interval B = 
12 — 34 T. The plotted energy is scaled as E/^/B to emphasize 
the Dirac fermion-like features in spectra (indicated by dashed 
vertical lines). Arrows denote transitions arising at the K 
point which evolve (nearly linearly) with B. Data above B = 
23 T were taken on second sample with a higher density of 
graphite flakes. 



behind the sample and cooled down to a temperature 
of 2 K, was analyzed by a Fourier transform spectrome- 
TBl 16| . The transmission spectra were normalized by 



ter 

the transmission of the tape and by the zero-field trans- 
mission, thus correcting for any magnetic field induced 
variations in the response of the bolometer. The missing 
parts of the transmission spectra, indicated by grey areas 
in Figs. [Hand [21 correspond to the spectral ranges where 
the tape is completely opaque. 

Prior to presenting our experimental results, we out- 
line a simple model of bulk graphite based on SWM 
model [1, G3. Whereas the standard SWM model has 
seven tight-binding parameters 70, ... ,75, we limit our- 
selves here to only the most important hopping integrals 
70 and 71. In other words, we consider only the parame- 
ters which are relevant for the nature of the band struc- 
ture in a graphene monolayer and bilayer. In graphene, 
the intra-layer coupling parameter 70 is directly related 



to the Fermi velocity, c = -\/3ao7o / (2h) , where the atomic 
distance is cto = 0.246 nm [23|, and in a graphene bilayer, 
the inter-layer coupling 71 gives an estimate for the mass 
of the charge carriers, m — 7i/(2c 2 ). 

Within our simplified approach, we calculate the band 
structure along the H — K — H line of the Brillouin zone, 
i.e. for —0.5 < k z < 0.5, which is essential for electri- 
cal and optical properties of bulk graphite. For the in- 
plane dispersion of charge carries, we find Hi! 24 1 
that it has, for a given momentum fc z , the form of a 
graphene bilayer with an effective coupling A71, where 
A = 2 cos(7rfc 



|2fl[24| 



In a magnetic field, we obtain the LL spectrum for each 
effective bilayer, i.e. for each momentum k z : 
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where ± sign labels the electron(+) and hole(-) levels. 
LLs related to the touching electronic bands are obtained 
for /i = — 1 and those related to bands split-off in en- 
ergy by an amount ±A7i are represented by /i = 1. 
The touching bands can be in the parabolic approxima- 
tion characterized by the mass m = A7i/(2c 2 ). The in- 
plane coupling 70 enters Eq. |T]) via the Fermi velocity c 
and directly influences the characteristic spacing of lev- 
els Ei = cy/2ehB. Note, that our approach is a special 
case of the model used by Koshino and Ando [21| . who 
in an analogous way calculated the spectrum of multi- 
layer Bernal-stacked graphene with an arbitrary number 
of layers in an external magnetic field. 

The joint density of states (initial and final states), es- 
sential in our magneto-optical experiments, has in the full 
as well as in our reduced SWM model singularities at two 
distinct points of the Brillouin zone, at the K {k z = 0) 
and H (k z = 0.5) points, where electrons and holes are lo- 
cated, respectively. Hence, the magneto-optical response 
of bulk graphite should be governed by transitions be- 
tween LLs defined by Eq. JTJ) for A — > 2 (K point) and 
A — > {H point). Notably, there is no singularity for 
A = I, which corresponds to a real graphene bilayer. 
Consequently, bulk graphite should, in magneto-optical 
experiments, behave as a combination of a graphene bi- 
layer with the effective coupling 271 and of a graphene 
monolayer, but with a twofold degeneracy _ 1 = £^ +1 1, 
in addition to the usual twofold spin and valley degen- 
eracies. The expected magneto-optical response of bulk 
graphite should therefore contain hole-related features 
whose energy evolves linearly with \/~B originating in the 
vicinity of the H point together with absorption lines 
whose energy evolves roughly linear with B correspond- 
ing to electrons at the K point. 

The magneto-transmission spectra taken for magnetic 
fields B — 12 — 34 T on a thin layer of bulk graphite 
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FIG. 2: (color online) (a): Positions of the 
absorption lines related to the H point as a 
function of \[B. The solid and dashed lines 
represent expected positions of absorption 
lines for c = 1.02 x 10 6 m/s (70 = 3.2 eV). 
(b): K point related absorption lines as 
a function of B. The solid lines show 
expected dipole allowed transitions in a 
graphene bilayer with an effective coupling 
271 calculated using Eq. (UJ) for 70 = 3.2 eV 
and 71 = 0.375 eV. Grey data points were 
taken on highly oriented pyrolytic graphite, 
which exhibit a behavior nearly identical to 
natural graphite [2^. The inset schemati- 
cally shows the observed inter-band transi- 
tions in the effective bilayer. 



are presented in Fig. [TJ The transmission is plotted as a 
function of energy divided by \/B~ to facilitate the identi- 
fication of spectral features originating from around the 
H and K points. Plotted in this way, the transitions 
denoted by Roman and Greek letters, do not shift for 
spectra recorded at different magnetic fields. Thus they 
scale linearly with y/B, see Fig. [2ji, and are related to the 
H point. The second set of lines marked by vertical ar- 
rows shift with magnetic field and actually follow a nearly 
linear dependence with B, as can be seen in Fig. ^Bp. 

The transitions following a y/~B dependence, corre- 
sponding to massless holes around the H point, have 
been thoroughly analyzed in our previous work 3, 22j |. 
Whereas, the absorption lines denoted by Roman let- 
ters have their direct counterpart in spectra of true 
graphene [H G3, H S HE Eg , the Greek lines are in 
principle dipole forbidden in a pure 2D system of Dirac 
fermions. Nevertheless, these transitions can be consis- 
tently explained with the same selection rule An = ±1, 
when the twofold degeneracy of LLs, _ 1 = e^+i v at 
the H point of bulk graphite is properly considered. The 
Fermi velocity is extracted to be c — (1.02 ± 0.02) x 
10 6 m/s, giving a rather precise measure of the in-plane 
hopping integral 70 = (3.20 ± 0.06) eV in bulk graphite, 
which is the only parameter needed to describe all ab- 
sorption lines originating at the H point. 

Here we focus on transitions denoted by arrows in 
Fig.[TJ whose energy evolves nearly linearly with B. Tak- 
ing into account the selection rule An = ±1 and using 
the in-plane coupling estimated above to 70 = 3.2 eV, we 
can interpret all the absorption lines in Fig. [2]d as dipole- 
allowed transitions in a graphene bilayer with an effective 
coupling 271 with 71 = (375 ± 10) meV. The deduced 
value for 71 is in a very good agreement with results ob- 



tained for a real graphene bilayer 28|, 2!| as well as on 
bulk graphite 3^, 31 1 . Our data follow well the theoreti- 
cal predictions up to the highest energies, in contrast to 
the the deviation reported for a real graphene bilayer [3] • 
Hence, the electrons in the vicinity of the K point can 
be described with a reasonable accuracy using the model 
of a graphene bilayer, but with coupling strength twice 
enhanced as compared to a true bilayer. The strength of 
the coupling 27! directly reflects the long-range Bernal 
stacking of graphite along the c axis. The twice enhanced 
coupling in the effective bilayer can be understood using 
the example of semiconductor superlattices, which are 
three-dimensional but nevertheless strongly anisotropic 
systems resembling in some aspects the band structure 
of bulk graphite. The energy difference, Asas, between 
the bonding and anti-bonding state in a symmetric dou- 
ble quantum well is simply half of the miniband width of 
the superlattice created from the same wells [32j . 

Our results also show that the parabolic approxima- 
tion which is widely used for the touching electronic 
bands in a bilayer and which directly leads to LLs whose 
energy evolves linearly with magnetic field: _ w 
±hw c yjn(n + 1) 0, [r| is a good approximation only 
in the vicinity of the charge neutrality point. As can 
be seen in Fig. [2j(b) , the small deviation from a linear 
dependence, predicted at higher energies by the simpli- 
fied SWM model for an effective bilayer, is reproduced 
in our data. The bilayer character of K point elec- 
trons also explains the non-linear evolution with B of 
the magneto-reflection spectra published recently 
Nevertheless, within the parabolic approximation, we ob- 
tain an effective mass of K point electrons in graphite of 
Tn = 71 /c 2 w 0.063mo in good agreement with other cy- 
clotron resonance experiments [33] . This mass is a factor 



4 



of two higher when compared to a true bilayer Q • 

While, only two tight-binding parameters 70 and 71 
are required to obtain a reasonable description of the 
magneto-optical response of the K point electrons over a 
wide range of energy and magnetic field, the influence of 
the remaining hopping integrals 72, . . . , 75, merits some 
consideration. In general, additional coupling parameters 
should lift the electron-hole symmetry of the bilayer and 
the trigonal warping (73) should lead to a mixing of LLs, 
which in turn can give rise to additional dipole-allowcd 
transitions. 

The electron-hole asymmetry at the K point should re- 
sult in an energy difference of the transitions L_„ — »L n _|_i 
and L_ n _! — >L n . Indeed, some evidence for this split- 
ting is visible in the spectra, see e.g. the transition n = 3 
above B ss 20 T in Fig. [1] Nevertheless, this effect is 
relatively weak, comparable to the width of the absorp- 
tion lines. The electron-hole asymmetry at the K point 
of graphite seems to be somewhat weaker in compari- 
son to the observed asymmetry in a true graphene bi- 
layer @,0,[33. 

The influence of the trigonal warping 73 on the 
magneto-optical response of the bilayer has been dis- 
cussed by Abergel and Fal'ko [Hj]. They conclude that 
this parameter becomes important only in the limit of 
low magnetic fields giving rise to a completely new set 
of dipole-allowed transitions. Similar reasoning [36| ex- 
plains the observation of numerous harmonics of the cy- 
clotron resonance of electrons in bulk graphite at low 
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Our simplified model also neglects the hopping inte- 
gral 72 which is directly responsible for the semi-metallic 
nature of bulk graphite. This parameter leads to a finite 
width of the doubly degenerate E3 and consequently, the 
fc z -dependent LLs n = and n = — 1 Assuming a 

finite and negative value of 72 [33j], the K point transi- 
tion L_i — >Lo cannot be observed at any magnetic field 
at low temperatures, as the Fermi level remains close to 
the middle of n = and n = — 1 Landau bands even in 
the quantum limit in graphite. 

To conclude, the magneto-optical response of bulk 
graphite can, over a wide range of energy and mag- 
netic field, be understood within a picture of an effective 
graphene monolayer and an effective bilayer with a cou- 
pling strength enhanced twice in comparison to a true 
graphene bilayer. This finding is in excellent agreement 
with predictions of a drastically reduced SWM model 
which retains only two tight- binding parameters, namely 
the in-plane and inter-plane coupling constants 70 and 
71. It should be noted, that as the validity of the model 
is limited in the vicinity of the Fermi level, it is not 
useful, for example, for the interpretation of magneto- 
transport experiments. Nevertheless, bulk graphite re- 
mains a material of choice to study magneto-optical phe- 
nomena in systems with both massless as well as massive 
Dirac fermions. 
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